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(54) Title: APPARATUS FOR THE MEASUREMENT OF GRAVITATIONAL FIELDS 
(57) Abstract 

Apparatus for measuring gravitational fields comprising a superconducting string 
(1) fixed at both ends and forming part of a closed superconducting loop inductively 
coupled to two driving solenoids (Ldi. La)- Displacement of the string in response 
to a gravitational field is sensed by two magnetic flux transformers each comprising 
a signal coil and two pick-up coils (Lpi, Lp2). Pairs of pick-up coils lie in two 
perpendicular planes providing two independent channels of measurements. The two 
arms of each flux transformer are balanced to convert only the amplitudes of the 
string's antisymmetric natural modes into an output voltage. The output voltage of 
each channel is used to produce a feed-back current distribution (I y i, Iy2) proximate 
and parallel to the string. By adjusting the feed-back current, the effective relaxation 
time and resonant frequency of the first antisymmetric mode of the string can be 
adjusted, while leaving the symmetric modes unchanged, thus increasing the apparatus' 
sensitivity to gravity gradients. 
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ap paratus measurement of gravitat i onal fields 

This invention relates to the measurement of 
gravitational fields, particularly to gravity 
gradiometry, and more particularly to a method for 
measuring absolutely off-diagonal components of the 
gravity gradient tensor. 



The gravity gradient tensor is a two-dimensional matrix 
of the second partial derivatives of a gravitational 
potential, V, with respect to the cartesian co- 
ordinates, x, y, z, of some arbitrary reference frame. 
15 It represents how the gravity vector itself in each of 
these directions varies along the axes . 

Accurate absolute measurements of the components of the 
gravity gradient tensor 1^ = aP^V (ij = x,y,z), taken at 
20 some local coordinate frame OXYZ are very important to 
progress in the fields of geological prospecting, 
mapping of the Earth's gravitational field, and space, 
sea and underwater navigation. 

25 A method of absolute measurement of gravity gradient 

tensor components was invented first by Baron Roland von 
Eotvos as early as 1890, utilising a torsion balance 
with proof masses hung at different heights from a 
horizontal beam suspended by a fine filament. The 

30 gravity gradients give rise to differential, forces being 
applied to the masses which result in a torque being 
exerted on the beam, and thus to angular deflection of 
the masses which can be detected with an appropriate 
sensor. A sensitivity of about 1 E (1 E = 1 E6tv6s = 10- 

35 s- 2 ) can be reached but measurement requires several 
hours at a single position due to the necessity to 
recalculate the gravity gradient components from at 
least 5 independent measurements of an angular 
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deflection each with a different azimuth angle. 

Practical devices, which have been built in accordance 
with this basic principle, are large in size and have 
5 low environmental noise immunity, thus requiring 

specially prepared conditions for measurements which 
excludes any possibility of using them on a moving 
carrier . 

10 A method for absolute measurement of gravity gradient 
tensor components which enhances the above method was 
invented by Forward in the middle of the sixties (see US 
patents 3,722,284 (Forward et al) and 3,769,840 
(Hansen) . The method comprises mounting both a dumbbell 

15 oscillator and a displacement sensor on a platform which 
is in uniform horizontal rotation with some frequency Q 
about the axis of the torsional filament. The dumbbell 
then moves in forced oscillation with double the 
rotational frequency, whilst many of the error sources 

20 and noise sources are modulated at the rotation 

frequency or not modulated (particularly 1/f noise) . 
The forced oscillation amplitude is at a maximum when 
the rotation frequency satisfies the resonance condition 
2Q=G> 0 / where co 0 is the angular resonant * frequency, and 

25 the oscillator quality factor Q tends to infinity. 

Unlike the non-rotating method, this method enables one 
to determine rapidly the quantities Tyy - and by 
separating the quadrature components of the response 
using synchronous detection with a reference signal of 

30 frequency 2Q. 

The same principles can be directly used, as proposed by 
Metzger (see US patent 3,564,921), if one replaces the 
dumbbell oscillator with two or more single 
35 accelerometers properly oriented on such a moving 

platform. There are no new features of principle in 
this solution to compare with the previous one except 
that the outputs of the pairs of accelerometers require 
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additional balancing. 

Devices have been built according to this method, but--;, 
they have met more problems than advantages, principally 
5 because of the need to maintain precisely uniform 

rotation and the small displacement measurement with 
respect to the rotating frame of reference. The devices 
have reached a maximum working accuracy of about a few 
tens of Eotvos for a one second measuring interval, and 

10 they are extremely sensitive to environmental 

vibrational noise due to their relatively low resonant 
frequencies. The technological problems arising in this 
case are so difficult to overcome that the existing 
developed designs of rotating gravity gradiometers are 

15 so far only at the stage of prototypes whose measurement 
accuracy is much lower than the limiting theoretical 
estimates . 



20 



In a paper by A. Nicolaidis and A. Taramopoulos (II 
Nuovo Cimento, Vol. 107B, N.ll, pages 1261-1266, 
November 1992) , the theoretical motion of a string with 
fixed ends under the influence of a plane monochromatic 
time-varying gravitational wave is discussed. According 
to this document, a string with fixed ends may be 
25 excited to resonance provided certain conditions, 

dependent on the length and orientation of the string 
and the wavelength of the gravitational wave, are met. 
It is suggested that Fourier analysis of the motion of 
the string could be used to extract the direction and 
30 energy of the incident wave. The document specifically 
avoids any discussion of the technical implementation of 
the theory, but it does suggest that strings a few 
metres or a few kilometres long should be used for the 
detection of cosmological radiation or gravitational 
radiation from a black hole or supernova, as the length 
of the string should be comparable to the wavelength of 
the gravitational waves. For the theoretical detector 
to work requires the gravitational field to oscillate in 



35 
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the form of a gravitational wave, which would not be the 
case for the gravitational' fields associated with 
massive bodies such as the Earth. 

5 Superconducting gravity gradiometers are known (see US 
patent 4,841,772 and Australian patent application 
48185/90) utilising a pair or more of sufficiently- 
separated superconducting accelerometers . Even after 
greatly reducing the intrinsic and environmental thermal 

10 noise factor, using stable persistent super- currents to 
balance the outputs of the accelerometers, and the most 
sensitive displacement sensors based on SQUIDs 
(Superconducting Quantum Interference Devices) , they 
cannot measure the gravity gradient tensor components in 

15 their absolute units because they are incapable of 
fixing a position of the accelerometer ■ s proof mass 
which is free of all forces. Therefore, only relative 
displacements of the proof masses can be measured. 
Rotating designs of such superconducting gravity 

20 gradiometers are not known. 

It is an object of the present invention to provide an 
apparatus for the measurement of gravitational fields 
with improved sensitivity, portability and noise 
25 immunity over the above known systems. 

It is a further object of the present invention to 
provide a novel apparatus for the absolute measurement 
of off -diagonal components of the gravity gradient 
30 tensor, in which the effect of rotation is replaced by 

parametric interaction between the sensitive element and 
active force feed-back connections, whereby enhanced 
sensitivity and vibrational noise immunity are attained. 

35 It is another object of the present invention to provide 
a simple technological realisation of the above 
apparatus utilising the advantages of the standard 
superconducting techniques which have shown an ability 
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to reach a maximum sensitivity for mechanical 
displacement measurements' and to keep intrinsic noise.at 
a minimum level- " 

5 to achieve these objects the present invention provides 
an apparatus for the measurement of quasi-static 
gravitational fields, comprising: 

a string fixed at both ends and held under tension; 
sensing means for detecting the transverse displacement 
10 of said string from an unperturbated position, due to a 
gravitational field acting on said string; 

and means responsive to the detected displacement 
for producing an output which is a function of the 
gravitational field. 

By -string" no particular limitation as to material or 
construction is intended. Any elongate tension element 
is included which is capable of being transversely 
deflected by a gravitational field and of providing a 
20 restoring force. 

An unperturbated stretched flexible string with fixed 
ends forms an absolute straight line in space going 
through the points where the ends of the string are 

25 fixed. This line can be identified as one of the axes of 
the local coordinate frame, say Z, and the other two 
axes, X and Y, are chosen to lie in the transverse to 
the string) plane. Any string deflection from thxs Ixne 
is caused by absolute values of the transverse 

30 components of the force per unit length which is applxed 
to each unit element of the string. 

Viewed from another aspect the invention provides a 
method of measuring quasi-static gravitatxonal fxelds, 

" C ° raP tro:iding a string with fixed ends and held under 

tension; m aa 
detecting the transverse displacement of saxd 
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string from an unperturbated position due to a 
gravitational field acting on said string; and 

producing an output in response to the detected , 
displacement, said output being a function of said 
gravitational field. 

The string -s deflection from its unperturbated position 
can be easily detected, by any suitable displacement 
sensing device. 



Preferably the string is formed of conductive, most 
preferably superconductive material. In this case, if 
an electric current flows through the string, a magnetic 
field distribution is produced in the transverse plane 
15 and along the string's direction. If the string is made 
of a superconducting material, a maximum current can be 
carried, and a consequent maximum sensitivity to the 
deflection can be reached. A d.c. or an a. c. current 
may be produced in the string by incorporating the 
20 string into a current -carrying circuit directly or by an 
inductive coupling with a pumping circuit (s), provided 
that the string forms part of a closed conducting or 
superconducting loop. An a.c. current may be induced in 
the string, for example by means of one or more, 
25 preferably longitudinally symmetrically positioned, 

coils, which may possibly be superconducting. The use 
of an a.c. current is advantageous in that it allows 
synchronous detection of the output signal. 

30 When the string carries a current, the transverse 

magnetic field around the string may interact with other 
conductors, or superconductors, by inductive coupling. 
The amplitude of the current induced in a conductor 
adjacent the string will be directly related to the 

35 distance of the string from that conductor. Thus, in a 
preferred embodiment of the invention one or more fixed 
pick-up coils are arranged along the length of the 
string to act as displacement sensing means, the current 
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induced in each coil being directly related to the 
string's displacement from its unperturbated position. 

In a preferred embodiment of the invention the sensing 
5 means comprises at least two sensors, possibly pick-up 
coils, positioned symmetrically, in the longitudinal 
direction, with respect to the mid point of the string. 

In a particularly advantageous embodiment, displacement 
10 sensors, for example pick-up coils, are arranged 

adjacent the string in two non-parallel preferably 
orthogonal, planes, so as to be capable of measuring the 
string's displacement in two transverse directions 
simultaneously . 



15 



20 



25 



30 



35 



It will be understood that the displacement of a string 
of length 1 from its unpertubated position, for example, 
in the y-direction. of the above local coordinate frame 
as a function of the z -position of a unit element and 
time, y(z, c) , can be described by the following 
differential equation 

--ng y (o,t)-nr yx (o,t)*./ i; (x, t> <D 

with boundary conditions corresponding to the fixed ends 
of the string, i.e. y(0,t) = yd, t) = 0. In this 
equation n denotes the string's mass per unit length, h 
is the friction coefficient per unit length, the 
parameters Y, A and Al/1 are the string's Young modulus, 
the area of its cross section and the string's strain 
respectively. The quantities g y (0,t) and r„(0,t> are the 
absolute values of the y-component of the total 
acceleration and the corresponding gravity gradient 
tensor component along the string, both taken at the 
centre of the local coordinate frame chosen. The 
function f>,t) represents the Langevin random force 
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25 



(2) 



per unit length acting on the string due to its 
interaction with the thermostat having the absolute 
temperature T, with the following correlation function^ 

5 ( f L iz,t)f\(z'.t')) - 4k B Th5 (z-z') 6 <t-fc'> 

where k B = 1.4 10" 23 JK" 1 is the Boltzmann constant and 
6(x-x') is the delta-function. 

10 in this description, the y-direction has been chosen as 
an arbitrary example to simplify the explanation of the 
invention. However, the foregoing and following 
analysis is equally applicable to any direction 
transverse to the string or any number of directions. 

15 

Applying Fourier analysis to the complex shape of the 
string caused by its interaction with the gravitational 
field, the function y(z,t), can be described, in the 
range z=0 to by an infinite sum of sinusoidal 

20 functions of period 21, with appropriate coefficients 

c y (n, t) . Thus a solution of Eq. (1) , which satisfies the 
boundary conditions shown above, can be represented by 
the following sum wherein each term in n corresponds to 
one of the string's natural vibrational modes 

y, Z ,t,.Ec y( n,t,si„(^) C3) 

By substituting Eq. (3) into Eq. (1) and by multiplying 
its left-hand and right-hand sides by sin (nn -z/1) , and 
then by integrating both sides over z from 0 to 1, one 
can obtain the differential equation for Cy(n,t) 



30 



•£c y (n,t>. 2^c y (n,t>^c y <n, t >- -L[<-1> » -l]g y <0, t>. 

35 21 

♦ (-1) "4^(0,0 



nn 



. -l-Jdz/ t U,t)sin^ip-zj ( 4 ) 
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20 



25 



30 



35 



where the quantities 



nn 

CO • A 

" 1 \ 



r ai (5) 

p 1 



represent the string's natural frequencies; x and p are 
the relaxation time and the volume mass density of the 
string respectively. 



When n takes an even value, i.e. for those terms of the 
infinite sum in Eq. 3 corresponding to vibrational modes 
of the string having a node at z=l/2 (antisymmetric 
modes), the term involving g y (0,t) is equal to zero. 
15 Thus, for n even, c y is dependent only on r yz (and 
thermal noise) . 

In practice this means that the amplitude, c y , of the 
antisymmetric sinusoidal components of the deflection of 
the string in the y-direction, yU, t) , is dependent only 
on the magnitude of the gravity gradient tensor 
component r^. 

The mid point of the string, z=l/2. is the position of a 
node in all antisymmetric vibrational modes of the 
string. If sensors are positioned symmetrically in the 
longitudinal direction with respect to this point, it 
will be possible to identify displacements of the string 
corresponding to the string's natural antisymmetric 
vibrational modes while discounting displacements 
corresponding to symmetric modes, the magnitude of which 
is not only affected by the gravity gradient tensor 
component T yx but also the absolute acceleration due to 
gravity in the y-direction, g y . 

It is particularly advantageous if displacement sensors 
are positioned at z=l/4 and z=3l/4, positions 
corresponding to the antinodes of the first 
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antisymmetric vibrational mode of the string, n=2. At 
these points the displacement of the string 
corresponding to the n=2 mode is at a maximum and thus-, 
the sensing signal will also be at a maximum, giving 
optimum sensitivity. 

According to a further development of the invention a 
conductor may be provided adjacent the conductive 
string. The conductor may carry a current directly 
related to the output of the sensing means, by the use 
of a positive feedback loop. The current may be 
activated continuously or periodically, for example in 
an -off -on" manner. In this case, a small deflection of 
the string due to a gravitational field will be 
15 amplified by the magnetic interaction of the string and 
conductor. In other words, the conductor will "push" 
(or "pull") the string into further deflection in direct 
response to a small deflection caused by a gravitational 
field acting on the string. This is clearly 
20 advantageous in that the displacement of the string is 
greater by virtue of the magnetic interaction with the 
conductor and is therefore more readily measurable, 
improving the sensitivity of the apparatus. 

25 In a particularly advantageous embodiment of this 
development, two or more conductors, possibly 
superconductors, are positioned longitudinally 
symmetrically about the mid-point of the string so that 
they amplify the antisymmetric modes of the string. 

30 

In overview, a preferred embodiment of the invention 
provides a novel apparatus for measuring absolutely and 
simultaneously a pair of off-diagonal components of the 
gravity gradient tensor by means of a flexible 
35 superconducting current -carrying string with fixed ends, 
comprising active parametric force feed-back 
connections. The string is the coherent sensitive 
element whose symmetric natural transverse modes are 
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caused by the total acceleration in the transverse 
plane, whilst the antisymmetric modes are caused only by 
absolute values of the gravity gradient components aldng 
the string's direction. 

in this embodiment the string forms a low- inductance 
part of a closed superconducting loop which is 
inductively coupled to a high- inductance driving 
solenoid(s) carrying an a.c. reference current from an 
10 external pumping source with some frequency Q. The 
string is also inductively coupled to two 
superconducting magnetic flux transformers each 
comprising a signal coil and two pick-up coils wherein 
the pairs of pick-up coils lie in two perpendicular 
15 planes the cross-line of which coincides with the 
unperturbated string, thus forming two independent 
channels of measurements. The two arms of each 
superconducting flux transformer are balanced to convert 
only the string's antisymmetric natural modes into 
20 signal current in the signal coil to be measured with 
SQUID • s (Superconducting Quantum Interference Devices) 
electronics. The output voltage of each channel, deeply 
modulated with the frequency O, is then proportional to 
the amplitudes of the antisymmetric natural modes of the 
25 string. This voltage is passed through a 

differentiating and summing amplifier, and then used to 
load the feed-back circuit to produce an in-line feed- 
back current distribution proximate and parallel to the 
string. By adjusting the feed-back current, the 
30 effective relaxation time and the resonant frequency of 
the first antisymmetric natural mode of the string 
(whose amplitude depends only upon the gravity gradient 
along the string's direction) can be increased and 
decreased respectively, while the same parameters, for 
35 the symmetric natural modes (whose amplitudes depend 

upon the total acceleration in the transverse- to-strmg 
plane) are not changed. In use, the feed-back circuit 
shifts the Brownian and vibrational noise level to far 
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below the sensitivity required for industrial 
applications . 

A preferred embodiment of the invention will now be 
5 described by way of example only and with reference to 
the following drawings in which: 

Fig. 1 is a general schematic representation of a 
preferred embodiment of the invention; and 

10 

Fig. 2 is a diagrammatic vertical cross-section of 
a device according to a preferred embodiment of the 
invention. 

15 A single channel prototype of a device according to the 
invention (see Fig. 2) has a flexible string 1. The 
string is preferably formed of a superconducting 
material such as Niobium (Nb) . Niobium wire is the best 
choice, having optimum elastic properties, which have 

20 been proven to be usable at 4.2 K. The string forms the 
low- inductance part L 0 of a superconducting closed loop 
which is inductively coupled to high- inductance driving 
solenoid (s) L d carrying an a.c. reference current I d (t) 
from an external pumping source with some frequency Q. 

25 The rest of the loop is provided by the casing of the 
device 2 , 2 1 , 3 , 4 , 5 . 

The string has, in this embodiment, a length 1 = 24 cm, 
is 1 mm in diameter and is fixed at its ends by two Nb 

30 cups 2,2' of cylindrical shape each having a hole of 1 
mm diameter at its centre. The cups 2,2' close tightly 
the ends of a Nb cylinder comprising three parts 3,4,5 
connected together with two Nb cylindrical rings 6,6» 
carrying a fine thread. The parts 3 and 5 also carry 

35 threads to engage other elements of the construction. 

The string's tension is produced by two Nb nuts 7,7' of 
1 mm fine thread. 
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The whole construction forms a closed superconducting 
cylindrical cavity with the string axially positioned. 
There are three spaces 10, 11. 12 inside this volume s 
electromagnetically insulated as much as possible from 
each other by Nb partitions 9. In two of them 10, 12, 
driving toroidal solenoids L dl and L^, wound with 0.01 ; 
mm Nb wire and connected in series, are placed, thus 
forming a large mutual inductance M d between. L^L^ + 
and the inductance of the cylindrical cavity L 0 which is 
of the order 10" 7 H for the sizes chosen. The ratio M d /L 0 
is about 5xl0 2 , so if the a.c. pumping current I d (t) in 
the driving solenoids has an amplitude of about 100 mA 
then the induced a.c. supercurrent I 0 carried by the 
string is about 50A peak to peak. In this case, the 
15 corresponding circular component of the magnetic 

induction B at the string's surface is nearly 200 Gauss 
which is approximately four times smaller than the 
Niobium first critical field B cl . 

20 The two rectangular- type pick-up coils L pl and L p2 of 
the superconducting flux transformer, and the active 
force feed-back circuit are mounted together on a 
Titanium tube 8 placed inside the central space of the 
construction shown in Fig. 2. Titanium is chosen because 

25 it has a thermal expansion coefficient matching that of 
Niobium. The active force feed-back circuit comprises 
two arms of 0.5 mm insulated copper wire stretched 
parallel to the string and carrying the feed-back 
current I y - I y i + lya* 

30 

This design has particular advantages; for instance, the 
closed superconducting configuration gives optimum 
shielding against external varying electromagnetic 
fields. Also, the cylindrically symmetric configuration 
35 has a small radial size which, including all integral 

parts of the prototype, is no more than 3.8 cm diameter. 
Thus, it is possible to utilise a standard commercial 
100 litre liquid helium vessel having an input opening 
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of about 4 cm diameter to cool the construction down 
with a standard probe. Special helium cryostats, which 
have been used for known devices, exclude the T % 
possibility of removing the device from the cryostat's 
5 inner volume, for example if something goes wrong, to 

readjust it under field conditions. Removal of a device • 
from a cryostat requires a long period, of for example 
several hours, to warm the cryostat contents to 
atmospheric temperature so that the contents do not 

10 explode under rapid thermal expansion. This is one of 
the major disadvantages of known constructions. 
However, the small input opening of the standard 
commercial 100 litre liquid helium vessel prevents such 
an explosion occurring, which means that the apparatus 

15 according to the present invention can be removed from 
the vessel and adjusted under field conditions. 

The string is deflected by a non-uniform quasi-static 
gravitational field and interacts with a variable feed- 

20 back current distributed close to and substantially 
parallel to the string. The distribution is optimum 
when the feed-back current is injected at or taken from 
the point in the feed-back circuit opposite the mid- 
point of the string (see Fig. 1) . In the local 

25 coordinate frame chosen this point is z = 1/2. Another 
requirement for optimum operation is that the two arms 
of the feed-back circuit are substantially equal and 
grounded at their ends. In this case, there is no 
electromagnetic coupling between the feed-back current 

30 and the closed superconducting loop in which the string 
is incorporated. 

The current I 0 (t) flowing through the string and 
interacting with the feed-back current distribution 
35 I y (z, t) gives rise to the following transverse component 
of the force per-unit -length f y (z, t) acting on the string 
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where Mo - 4n 10" 7 Hm" 1 is the magnetic vacuum 
5 permeability, d is the distance between the centre of 
the unperturbated string and the centre of the wire 
carrying the feed-back current and the phase of the 
pumping current source 15 is chosen to be zero. The 
sign + or - is determined by the output buffer of the 
10 differentiating and summing amplifier 16 shown in Fig.l. 
The transverse motion of the unit element of the string 
in the OYZ plane is then described by the following 
differential equation 



at 



.- n<r y <0,t)- nr^O, t)z i^V,l t) sin (Qt)- f t (z,t) 



(7) 



which will be seen to correspond to Eq. (1) with the 
addition of the term of Eq.(6) . Consequently, Eq. (7) 
also has solutions of the form of Eq. (3) . Thus, 
following the same algebraic manipulation as for Eq. (1) , 
a differential equation for c y (n, t) can be obtained for 
this embodiment. This is Eq. (8) which corresponds to 
Eq. (4) but with the addition of a feedback term. 



i-i) n — r (o, t) 

nn rt 



± lJt2-e n I 0 I v [z, t) sin (Qt) 
il 2n 2 d 

« JLfdzf L [z t t)sin^ij-zj (8) 

The quantity e n relates to the characteristics of the 
transducer system of the feedback loop; the longer the 
length of the arms of the feedback circuit, the larger 
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the quantities e n are. If the arms of the feed-back 
circuit are absolutely identical then the quantities e n 
are equal to zero for all odd n = 1,3,5.... Their 
particular values for the sizes shown in FIG. 2 are 
5 determined by 

So, for the properly adjusted feed-back circuit, only 
10 the antisymmetric natural modes of the string interact 
with the feed-back current. However, only the 
antisymmetric natural modes of the string are sensitive 
to the absolute value of the gravity gradient tensor 
component to be measured, as is seen from Eq. (4) and 
15 (7). 

The superconducting pick-up coils L pI and L p2 are placed 
near the string and cause two arms of the 
superconducting magnetic flux transformer to convert, if 

20 perfectly balanced, only the antisymmetric natural modes 
into the signal current Ii 'to be detected with the 
SQUID'S electronics 13 (see FIG. 1) . One uses SQUID'S 
(Superconducting Quantum Interference Devices) 13 as 
they are the most sensitive variable current and 

25 magnetic flux sensors currently available. In the 

prototype shown in FIG. 2, the pick-up coils are made in 
the form of two rectangular- type single loops of Nb wire 
placed symmetrically with respect to the midpoint of the 
string and connected in parallel with the signal coil Z^. 

30 If the symmetry is perfect and the areas of the loops 
are absolutely identical, then the symmetric natural 
modes do not produce any signal current I ± or feed-back 
current J y . The same effect can be achieved for slightly 
non- identical pick-up coils with the accuracy required, 

35 if one uses the additional inductance (s) L b connected in 
parallel and/or series with one or both of the pick-up 
coils. The inductance (s) L b can be. tuned to balance the 
two arms of the superconducting flux transformer. The 
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residual "zero-mode" current in the signal coil L £ 
corresponding to the unpefturbated position of the 
string can be compensated directly inside the SQUID bx 
an additional coupling (not shown) to the pumping 
current source. If the balancing conditions are 
satisfied, then the output voltage of the SQUID'S 
electronics 13 is determined by 

V y {t). KI 0 L a sin(Qt)|Ep n c y (n,t)j» K* M (t) (10) 



where K is the total flux to voltage transfer function 
and L s is the SQUID'S inductance. The quantities p n 
depend on the physical design and position of the pick- 
up coils and are equal to zero if n = 1,3,5 The 

15 function ^(t) is the equivalent-to-noise random magnetic 
flux inside the SQUID loop, whose spectral density S*(g>) 
determines the intrinsic instrumental limit of the 
accuracy of measurements. The feed-back current I y (t) is 
formed from the output voltage V y {t) by passing it 
through a differentiating and summing amplifier 16 which 
is loaded by resistance R y . In this case, the feed-back 
current I y (t) can be represented by 



where p, q and x* are constant parameters which depend 
upon the design of the differentiating and summing 
amplifier 16. 

30 It must be noted that a mismatch between the two arms of 
the feed-back circuit always exists. The design shown 
in FIG. 2 uses two identical feed-back resistances, R yl 
and R y2 , one for each arm. In this case the mismatching 
can be easily compensated by tuning one of the 

35 resistances, say R y2 , to obtain the optimum case. 

Equations (7) and (11) represent a closed infinite set 
of differential parametric- type equations. Careful 
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analysis has shown that we can ignore the terms 
involving the quantities' c y (n, t) with n> 2 in the right- 
hand side of Eq. (11) . The reason is that just one mqde 
can be made "soft" i.e. the most sensitive to the 
gravity gradient, namely c y (2,t). If the string's 
natural frequencies are high enough and separated by 
single octave gaps, only second-order corrections are 
required which can be easily taken into account along 
with analysis of other instrumental errors. Then, as 
follows from Eq. (7) , the self -consistent equation for 
the gravity gradient sensitive mode, n=2, including 
unavoidable fundamental noise sources can in turn be 
written in the form 



d 2 

15 ~T~i c y 



( 2,t)^l - {a - jacos(20t))-jU y <2,t). 
.|c)2_ 1<5K io5 2 cos (20t)- i.aQsin(2Qt)jc y (2,t)- JLr„{0,i). 

, JLfdzf (z, t) sinf — z\* (back-action noise ) sin (Qt) (12) 

where 

(13) 



25 and it is assumed that the . true sign of the feed-back 
current has been chosen. 

If some easily carried out conditions are satisfied, 
which are 



1 dT 
r dt 



a 2. |S 2 8Q 2 »co 2 < 14) 



then one can show that the self -consistent output 
voltage is 



35 



v (t)-B,KX 0 i 



f 1 ^(Ort) 

1 11 * T S2 



brownian noise 



sin (Qt)- K* N [ t) (15) 
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where under "brovmian noise" the combination of thermal 
and back-action noises is implied. 

■d& ■ 

It is of interest to estimate the limiting accuracy of 
measurements of this embodiment of the invention, which 
can be represented by the value of a minimum detectable 
gravity gradient 



10 IW 



ii 
1 



16Jc fl T J5< VQ> ll?9 Eotvds . ( 16) 

4 n\ ^ 



where T eff = x/(l-cn:/4) is the effective relaxation time, 
m is the total mass of the string, and £*(Q) is the 
15 energy resolution of the SQUID. Using the following 

practical parameters: 1=0. 24m, m-1.6 lO' 3 kg, x eff - 10 4 s, 
p 2 - 4xl0 3 m- 1 , L s -5xl0-" H, I„- 50 A, (a) 2 2 - U 2 /2) 1/2 /2n« 2 Hz, 
co 2 /2n- 40 Hz, Q/2H i 2xl0 2 Hz, E^Q) - 2xl0" 31 J/Hz (d.c. 
biased SQUID), one can obtain from. Eq. (16) 



n * Eotvos (17) 
r »m" °- 4 — 7=^- 

It can be shown, that a range of the parameters x, o> 2 , ca, 
a and D exists where the string's response described by 
25 Eq. 12 is stable. For example, for quasistatic gravity 
gradients and sufficiently high pumping frequency O one 
can ignore the oscillating terms containing Cos(2Qt) and 
Sin(2Qt) in the right side of Eq. 12. 



There are a number of detecting strategies which can be 
employed by the present invention at this stage, which 
are dependent on the initial mechanical parameters of 
the string and the application for which the apparatus 
is intended. It is preferable to use a string with a 
high mechanical stiffness and a short relaxation time in 
order to increase immunity to vibrational noise, which 
is the main noise source in industrial applications, 
particularly in mobile gravity gradiometry. On the 
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other hand, the stiffer the string, the stronger the 
feedback force that has to be applied to the string to 
soften the signal mode, and the larger the back-actiSh 
noise associated with the feedback current. 

Additionally, the shorter the string's relaxation time, ■ 
the stronger the influence of thermal fluctuations of 
the string on the measuring accuracy since the mass per 
unit length of the string will normally be quite small. 

To overcome both of these problems a best mode of 
carrying out gravity gradient measurements according to 
another embodiment of the present invention uses 
variable feedbacks in an "off-on" manner. In this case, 
the feedback force is initially not applied to the 
string for an 'off -period' during which the string 
reaches thermodynamic equilibrium. The feedback force 
is then quickly activated for an 'on-period' during 
which the effective natural frequency 



and the effective relaxation time 



become substantially smaller and longer respectively 
compared to the corresponding initial parameters of the 
string. The feedback is adjusted in such a way that the 

30 effective relaxation time becomes much longer than the 
on-period. Measurements are carried out during the on- 
period only, in which the string never reaches 
thermodynamic equilibrium. For example, the fluctuation 
dissipation theorem is no longer applicable to the 

35 string during the period of measurements and its 

response to all external noise sources is changed (see 
V.B. Braginsky and A.B. Manukin, Measurement of Weak 
Forces in Physics Experiments, Ed. by D.H. Douglass, 
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University Press of Chicago, 1977) . 



One can show that in this case the least gravity 
gradient detectable by this embodiment of the invention 
can be represented by 



m n 

mln " 2 



where 



16 V *» (Q) 1h( —I 10» Eotvos (20) 



v. 6 .-^m)-^[i-(i)] 



(21) 



T m is the measurement time (on-period) , m is the total 
mass of the string, E*{Q) is the energy resolution of the 
SQUID at the frequency Q and 6 is a statistical error of 
the first kind. The value of 6 is the likelihood that 
the equivalent gravity gradient noise will exceed the 
value represented by the left side of Eq. 20 for the 
period of measurement. 

Using the following practical parameters: 1=0.24 m, 
ro =1.6xl0- 3 kg, x=0.5 s, x m =l s, x e££ « 10* s, p 2 - 4xl0 3 m" 1 , 
L s "5xl0 _11 H, J 0 - 50 A, G> ef£ /2n« 3 Hz, o 2 /2n« 80 Hz, Q/2n s= 
10 4 Hz, MQ)- 5xl0- 32 J/Hz (5001* d.c. biased SQUID), one 
can obtain from Eq. (20) 



30 =0.02 Eotvos 



In both the above embodiments, the desired signal is 
obtained from the output voltage by synchronous 
detection with a reference signal taken from the pumping 
35 source 15, and the invention allows calibration of the 
desired signal in gravity gradient absolute units 
without rotation as has been proposed for known rotating 
gravity gradiometers . As for rotating designs, the 
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invention allows the movement of the noise spectrum to a 
frequency range at which' 1/f contribution is 
sufficiently small* Natural vibrations of the string, 
which occur during the time of measurement (on-period) , 
5 do not cause a problem since they can be filtered out 
from the desired signal provided that the on-period is 
chosen to be much longer than the period (2n/G> eff ) of 
such vibrations. 

10 Vibrational noise immunity is improved by the factor 
(cOeff/G)!) 2 which can be made as small as 10~ 2 . 

One must consider inductive cross-coupling between the 
feedback currents and each pair of the pick-up coils and 

15 cross -coupling between the pick-up coils themselves, 

both of which act like negative feedback loops. On the 
one hand this leads to unnecessary renormalisation of 
the amplitudes of the output signals until the gain of 
the SQUID'S electronics exceeds some critical value. On 

20 the other hand in the case of double-channel 

measurements, the output signal of each channel contains 
a linear combination of each gravity gradient component 
to be measured. It can be shown that each of such 
components can, nevertheless, be measured separately and 

25 simultaneously, if a proper data acquisition system is 

used. The effect can be easily eliminated by organising 
additional positive feedback to counteract this negative 
feedback, for example by connecting, via a weak 
inductive coupling, each feedback current with each 

3 0 SQUID. 

In practice, the apparatus according to the invention 
can be used to determine in absolute units the off- 
diagonal components of the gravity gradient. By 
35 conducting a gravity survey over an area, small 
differences in absolute gravity gradient can be 
detected. Such small changes may indicate variations in 
local geological features, for example, the presence of 
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minerals, gas or oil. 

Repeated readings over time at a single locality could 
indicate changing geological status of an area, such as 
rising magma. Clearly the invention enhances 
prospecting and other data gathering pursuits where 
accurate gravitational field measurement is required. 
Use of absolute values enhances the information that can 
be determined from the data measured. A gradiometer 
according to the invention can be used while moving, 
which allows the gradiometer to be used on vehicles 
whether land, sea or air vehicles. For example, the 
device can be suspended from a helicopter and used while 
the helicopter traverses a selected area. 
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Claims 

1. Apparatus for the measurement of quasi-static „ 
gravitational fields, comprising: 

5 a string fixed at both ends and held tinder tension; 

sensing means for detecting the transverse 
displacement of said string from an unperturbated 
position due to a gravitational field acting on said 
string; and 

0 means responsive to the detected displacement for 

producing an output which is a function of the 
gravitational field. 

2. Apparatus as claimed in claim 1, wherein said 
5 sensing means comprises at least two sensors 

symmetrically longitudinally positioned about the mid- 
point of said string. 

3. Apparatus as claimed in claim 1 or 2, wherein said 
0 string is composed of conductive material and carries a 

current J 0 . 

4. Apparatus as claimed in claim 3 further comprising 
conductive means adjacent said string carrying a current 

5 I yt wherein the magnitude of the current I y is a function 
of the output of the sensing means; and the magnetic 
fields associated with the current I y through the 
conductive means and the current J Q through the string 
interact to produce a feedback force on said string so 

0 as to increase the transverse displacement of the string 
from its unperturbated position in response to the 
gravitational field acting on said string. 

5. Apparatus as claimed in claim 4, wherein said 
5 conductive means comprises at least two conductors 

longitudinally symmetrically positioned about the mid- 
point of said string, each conductor carrying a 
substantially equal proportion of said current I y . 
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6. Apparatus as claimed in claim 4 or 5 wherein the 
current I y through said conductive means is activated 
periodically. 

7 . Apparatus as claimed in any of claims 3 to 6, 
wherein said sensing means comprises at least one pick-, 
up coil in which a current J p is induced by said current 
J c through said string, said current J p being a function 
of said string's displacement. 

8. Apparatus as claimed in any of claims 3 to 7, 
wherein said current I 0 through said string is an 
alternating current. 

15 9. Apparatus as claimed in any of claims 3 to 8 

wherein said current I 0 is induced in said string by 
inductive means. 

10. Apparatus as claimed in claim 9, wherein said 

20 inductive means comprises two solenoids longitudinally 
symmetrically positioned about the mid-point of said 
string . 

11. Apparatus as claimed in any of claims 3 to 10, 
25 wherein said string is composed of superconducting 

material . 

12. Apparatus as claimed in any preceding claim, 
wherein said sensing means comprises means for detecting 

30 the transverse displacement of said string in two non- 
parallel planes. 

13 . Apparatus for measuring of f -diagonal components of 
the gravity gradient tensor comprising a flexible string 

35 having fixed ends, means to apply a force per unit 

length to said string so that deflection of said string 
is caused by absolute values of transverse-to-string 
components of said force, in a manner such that ■ 
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deflection of said string is a combination of said 
string natural modes, and the even said modes are caused 
only by absolute values of the components of the gravity 
gradient in the string direction, whilst the odd said 
5 modes are caused by total acceleration in the 
transverse- to- string plane. 

14. A method of measuring quasi-static gravitational 
fields, comprising: 

10 providing a string with fixed ends and held under 

tension; 

detecting the transverse displacement of said 
string from an unperturbated position due to a 
gravitational field acting on said string; and 
15 producing an output in response to the detected 

displacement, said output being a function of said 
gravitational field. 

15. A method as claimed in claim 14 wherein said output 
20 is produced by measuring the spatial position of at 

least one point on the string relative to the 
unperturbated position of said point. 

16 . A method as claimed in claim 15 wherein the spatial • 
25 positions of an even plurality of points on the string 

are measured relative to their unperturbated positions 
and said points are chosen to lie pair-wise 
longitudinally symmetrically about the mid-point of the 
string. 

30 

17. A method as claimed in claim 16 wherein said points 
correspond to the positions of antinodes of the 
antisymmetric natural modes of said string. 

35 18. A method as claimed in any of claims 14 to 17 

wherein said displacement of said string is increased by 
applying a feedback force to the string, said force 
being a function of the gravitational field acting on 
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the string. 



19. A method as claimed in claim 18 wherein said 
feedback force is a direct function of said output. 

20 . A method as claimed in claim 19 wherein said 
feedback force is applied to said string so as to 
accentuate the components of the spatial configuration 
of the string corresponding to natural antisymmetric 

10 modes of the string in preference to the components of 
the spatial configuration corresponding to natural 
symmetric modes. 

21. A method as claimed in any one of claims 14 to 19 
15 wherein said displacement is measured in two non- 
parallel planes. 

22 A method of measuring absolute off -diagonal 
components of the gravity gradient tensor by means of a 

20 flexible string with fixed ends whose deflection from 
its unperturbated position is caused by absolute value 
of transverse-to-string components of a force per-unit- 
length applied to each unit element of said string m a 
manner that said deflection is bt combination of said 

25 string's natural modes, and the even said modes are 

caused only by absolute values of the components of the 
gravity gradient in the string -s direction, whilst the 
odd said modes are caused by total acceleration in the 
transverse-to-string plane. 



30 
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